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The Optimized Perturbation Theory (OPT) method, at finite temperature and finite chemical 
potential, is applied to the field theory model for polyacetylene. The critical dopant concentration 
in trans-polyacetylene is evaluated and compared with the available experimental data and with 
previous calculations. The results obtained within the OPT go beyond the standard mean field 
(or large-iV) approximation (MFA) by explicitly including finite N effects. A critical analysis of 
the possible theoretical prescriptions to implement and interpret these corrections to the mean field 
results, given the available data, is given. For typical temperatures probed in the laboratory, our 
results show that the critical dopant concentration is only weakly affected by thermal effects. 

PACS numbers: 64.60.A-,71.30.+h,11.10.Kk 



About three decades ago a remarkable discovery was made that trans-polyacetylene (GH) a doped with halogens 
could behave as a metal, exhibiting electrical conductivity [1]. Since then several striking features have been shown 
by conjugated polymers, such as electronic, optical and magnetic properties, which give these materials a wide range 
of applicability 0, 0, 0, Before reaching the metallic state, polyacetylene can be converted into a semiconductor, 
depending on the concentration of dopant, y, defined as the number of doped electrons per carbon atom. For lower 
dopant concentration, the conduction can be described in terms of topological excitations, such as (spinless charged) 
polarons, bipolarons or solitons [6(. Experiments have also shown that the observed non-metal to metal transition 
in polyacetylene, that typically happens when the dopant concentration is increased to a critical value, y c ~ 6%, 
corresponds to a first-order transition Q. 

From the theoretical point of view, the role of the (self-localized) solitons in the charge-transport process is suc- 
cessfully described by the model Hamiltonian proposed by Su, Shrieffer, and Heeger (SSH) Q. SSH were primarily 
interested in the low-energy excitations of the one-dimensional polymer. In this case, the electronic correlation 
length, £, becomes much larger than the lattice constant, a, and an effective model Hamiltonian incorporating only 
the elcctron-phonon coupling as an interaction term, and neglecting electron-electron Coulomb repulsion, was shown 
to be able to give a very good description of the system. To take into account electron-electron interactions, extended 
versions of the SSH Hamiltonian, such as Hubbard models [5|, [9(, have been later employed. 

The continuum limit of the SSH Hamiltonian is known as the Takayama-Lin-Liu-Maki (TLM) model [ic|, which is 
a field theory with two-flavor Dirac fermions. The TLM model gives a very good approximation to the discrete SSH 
model provided two conditions hold: the band gap (denoted as E gap ) between a electron bonding and anti-bonding 
states is much smaller than the 7r-band width (denoted as W) and the typical distance between the neighboring 
single-electron levels near the Fermi energy is much smaller than the gap. The first condition is automatically fulfilled 
for trans-polyacetylene, where from the typical measured values E gap /W ~ (1.8eV)/(10eV) = 0.18 [f|, while the 
second one is met provided the application is for a chain with a large size (Tlj . At the same time, in many conducting 
polymers the characteristic correlation length, £, is much larger than the lattice constant, a. These two quantities are 
related to E gap and W by £/a = W/ E gap j(| . 

Some authors have recognized that the corresponding Lagrangian of the TLM model is analogous to a field 
theoretical model for fermions with quartic self-interactions in 1+1 dimensions, also known as the Gross-Neveu (GN) 
model 13]. It is worth recalling that the large- N, or mean-field approximation (MFA), predicts that, at T = 0, the GN 
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model in 1+1 dimensions suffers a first-order phase transition at the critical chemical potential value, p c = Aq/v2, 
where A is the order parameter value at T = = p, UM . Other works making use of the GN model applied to the 
study of the properties of polyacetylene include Refs. [la. Il6j. However, to our knowledge, the issue of investigating 
the phase diagram and temperature effects in the dopant concentration in polyacetylene have not been addressed, in 
the context of the GN model, so far. 

Since we can relate the dopant concentration in terms of the chemical potential, the GN model works as an analog 
model that can describe, in a field theory language, the non-metal to metal phase transition observed in the trans- 
(CH) X . An early study of the transition from the (chiral symm etry broken) low-density to the (symmetry restored) 
high-density phase, has been performed by the authors of Ref. [17| using the GN model, by assuming fermions with 
N = 2 flavors, but considering only the leading order in a 1/N expansion. Despite this approximation, it was found an 
excellent agreement, of about 6%, between the theoretical and experimental critical doping concentration, y c . Later, 
the same authors have used a mixture of the thermodynamic Bethe ansatz (TBA) and the 1/N expansion at the 
next to the leading order (NLO) to evaluate finite N corrections to p c ,18]. Surprisingly, their NLO does not agree 
so well with the experimental results, deviating by about 20%. One of our main motivations in the present work is 
to shed some light into this rather peculiar behavior. After all, one expects that finite N corrections should improve 
convergence, especially in a case where N = 2, like it is the case in the theoretical model description of polyacetylene. 
With this aim we start by remarking that the critical dopant concentration, y c , is related to the critical density, p c , 
and the equilibrium space, a, between the x coordinates of successive CH groups in the undimerized structure by 
He = Pcd- However, in Ref. [l8j . the authors did not compute the finite N corrections to the critical density p c 
but rather to the value of the chemical potential p c . Nevertheless, they concluded that, for N — 2, p c (N) should 
change by about 20% since it is proportional to p c (N), which changes by the same amount. We believe that this 
conclusion is not necessarily right since additional N and/or coupling (A) dependent factors, H(\,N), could appear 
when obtaining p by deriving the pressure with respect to p so that p c {N) could be proportional to H(X,N) p c (N) 
where, for example, TL{\, N — > oo) — > 1 in the mean-field approximation. Second, p depends on parameter values 
such as Sgan and W that sets £/a = W/E gap but only the value £/a ~ 7, which corresponds to the original SSH 
parameters Q E gap = 1.4 eV and W — 10 eV, was considered in Refs. [13) Hi- We must recall that the original SSH 
parameter values, especially the one regarding E gapi have been further updated already in the Ref. [6|. 

Here, we shall investigate the same problem considering again the 1+1 dimensional GN model, but not restricting 
the analysis to the large- N approximation. On general grounds, one may eventually wonder on the relevance of going 
beyond the large N (or equivalently mean field) approximation in the present context, given that the equivalence 
between the original TLM model and the GN field theory was established strictly speaking at the mean field level: 
indeed, there are important phenomenological aspects and properties of the polyacetylene that are anyway not cor- 
rectly described by the sole TLM (for a review on both the success and limitations of the TLM/GN equivalence see 
e.g. Nevertheless, it appears sensible here to push the correspondence further in this direction at least because 

the TLM model is effectively a two-flavor only model: since some results and techniques to go beyond large N are 
available, these can be expected to improve the comparison of the GN model to polyacetylene data, as we shall see 
later. In this work we use the method of the optimized perturbation theory (OPT) [20], which is a nonperturbative 
method allowing to consider corrections going beyond the large-A^ (or mean-field) approximation. Our aim is to 
obtain relations for p c which include finite N corrections. For this we will reconsider Landau's free energy density for 
the 1+1 dimensional GN model that has been recently derived by some of the present authors in Ref. [2l|, where the 
OPT results for the phase diagram show how this method correctly improves over mean-field results, in accordance 
with Landau's theorem for phase transitions in one space dimensions. In the same reference one finds N dependent 
analytical equations for the order parameter and chemical potential which turn out to be crucial for the present 
application. A detailed description of how the method works can be found in Ref. [2l| and references therein. 

In the present application, we will indeed show that within the OPT one can obtain the critical dopant concentration, 
y c , in a consistent way, producing results which are in very good agreement with the experimental values for the trans- 
(CH) X when up-to-date parameter values are considered. For completeness however, we will examine several possible 
theoretical prescriptions to implement and interpret these corrections to the mean field results. More precisely, as we 
shall see, the corrections to the mean-field may be interpreted in different ways within the context of polyacetylene. 
This is due to the fact that its description in terms of the Gross-Neveu field theory model is strictly speaking only an 
approximation. We shall also pay attention to the non negligible experimental error on some of the relevant available 
data. Apart from considering finite N corrections to y c , we also estimate how thermal effects affect this quantity. 

It is worth mentioning that an extension to the 2 + 1 dimensional GN theory is straightforward and may be of use 
in the investigation of superconducting electrons in quasi- two-dimensional systems [22j . In connection with the 2+1 
dimensional GN model, the OPT has allowed [23| to redraw its phase diagram leading to the precise location of a 
tricritical point in the T — p plane and the discovery of a mixed "liquid-gas" phase that remained undetermined for 
almost twenty years, since the first mean-field results appeared |24l|. The compatibility of the OPT and Landau's 
expansion for the free energy density has been shown in Ref. [25]. In the condensed matter domain, higher order 
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OPT results have considerably improved early applications of the same method [26| producing some of the most 
precise analytical values for the shift of the critical temperature, AT C , of an interacting homogeneous Bose gas when 
compared to the ideal gas [13] ■ Finally, the OPT convergence has also been proved in connection with critical theories 

As far as thermal effects in the GN model at one space dimension and the application of these results to polyacetylene 
are concerned, a few comments are appropriate here. We recall that due to a well known no-go theorem [30l l3ll l32|. 
for a one-dimensional system at any finite temperature we should expect no phase transition related to a discrete 
symmetry breaking (in this case a discrete chiral symmetry in the massless GN model considered in this work). This 
is due to kink-like inhomogeneous configurations [33| that come to dominate the action functional, instead of just 
homogeneous, constant field configurations. This is to be contrasted to the phase transition observed since long ago 
in the GN model in one space dimension in the mean- field, large- TV approximation [Til ]. This result is explained by 
the way the thermodynamic and the mean-field approximation are performed. If the thermodynamic limit is taken 
before the mean-field approximation, those large nonhomogencous fluctuations dominate and the theorem is observed. 
However, if the mean-field approximation is considered first, the fluctuations are suppressed, thus seem to evade the 
no-go theorem. Since we are here applying the GN model as an effective analog model for the polyacetylene and this 
is in practice a finite size system, we do not expect the theorem to be completely observed here. In fact, a phase 
transition at finite temperatures is indeed observed and measured in the laboratory. Nevertheless, polyacetylene is 
a well-known system exhibiting a rich spectrum of nontrivial fluctuations, from solitons to polaron excitations [||. 
Therefore, we may expect not only homogeneous like configurations (like in the mean-field approach), but also that 
the inclusion of these excitations in any theoretical calculation in this model should be considered. In this context, 
for example in the GN field theory model, by accounting for kink-like configurations in the large- -/V approximation, 
the authors of Ref. [34[ found evidence for a crystal phase that shows up in the extreme T ~ and large fi part of the 
phase diagram, while the other extreme of the phase diagram, for large T and (i ~ 0, seemed to remain identical to 
the usual large- N results for the critical temperature and tricritical points, which are well known results [Tij for the 
GN model. In this work we will only consider homogenous vacuum backgrounds in our thermodynamical calculations 
applied to the polyacetylene. By comparing our results with the experimental ones we can roughly estimate the 
importance of soliton-like excitations in the system. From our results, we estimate that these effects are expected to 
be small in the context of applying the GN model as an effective analog model to describe the thermodynamics of 
polyacetylene at low (laboratory) temperatures. 

This work is organized as follows. In Sec. II we briefly present the TLM model and its relation to a four-fermion 
theory, which can be identified as the GN model. In Sec. Ill we review the computation of the free energy for the 
GN model by using the OPT method. In the same section, the temperature dependent density is obtained. The gap 
equation is used to set up the parameter values in Sec. IV. In Sec. V we show our phase diagrams for the trans- 
polyacetylene (CH) X both in the T — (i and T — p planes. The critical dopant density, at zero and finite temperatures, 
is considered in Sec. VI. Our concluding remarks are given in Sec. VII. 

II. THE TAKAYAMA— LIN-LIU— MAKI AND GROSS-NEVEU MODELS 

The Takayama-Lin-Liu-Maki (TLM) Hamiltonian is the continuum version for the original SSH model and it is 
given in terms of a fermionic field, ip, and a scalar field, A, representing the coupling of the electron gas to the local 
value of the dimerization and it is expressed by the Hamiltonian [To| 

Htlm = - r \ / dxA 2 (x) + V / dx^(x) [-ikv F a 3 d x + aiA(x)} %j)(x) , (2.1) 

ZirhVF Atlm J s J 

where the sum is over the spin states, Oi are the Pauli matrices, vf — kp/m is the Fermi velocity and Atlm is a 
dimensionless coupling defined by 



AtlM " ' (2 - 2) 

where a is the 7r-electron-phonon coupling constant of the original SSH Hamiltonian, K is the elastic chain deformation 
constant and to is the hopping parameter, which is expressed in terms of the Fermi velocity and the equilibrium space 
a between the x coordinates of successive CH groups in the undimcrized structure as t$ = Hvf /(% a). 

Note from Eq. (|2.1[) that a nonvanishing (constant) value for A leads to a mass term for the fermions, thus breaking 
the chiral symmetry exhibited by .Htlm an d opening an electronic energy gap in the system. The presence of a gap 
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prevents electrons to move to the conduction band and, thus, the system effectively behaves as a non-metal. The 
effect of the addiction of dopants to the system is to decrease the electronic energy gap, till it vanishes at some 
critical dopant concentration and the system starts to behave as a metal. In general, a kinetic term for the scalar field 
emerges when taking the continuum limit of the SSH model. However, we consider the usual adiabatic approximation 
of neglecting the lattice vibrations, valid for energies for the optical- phonons (given by Tiujq) smaller than the gap 
magnitude (2A). In particular, for typical values found for polyacetylene [6(, 2A ss 1.8 eV and hujQ « 0.12 eV, this is 
regarded as a valid approximation. 

The model described by Eq. (|2.1|) can easily be shown to correspond to a four-Fermi model if we eliminate the 
scalar field A from Eq. (|2.1[) . e.g. by using its equation of motion. Then, putting the TLM model in the Lagrangian 
density form one obtains 

£tlm = -» ; \ A 2 + ^ (ihd t - ihv F j 5 d x - 7 A) ip , (2.3) 

ZttTlvf Atlm 

where we have identified 75 = —03 and 70 = o\. Now, eliminating A from Eq. (|2.3[) upon using 71 = 102, as well as 
the usual relations between the Dirac matrices, leads to 

£tlm = V> (ihjodt - ihv F 7id x ) ip + "T^T^-P , (2-4) 

which is just a four-Fermi Lagrangian density corresponding to the massless GN model [l3| where N is the number 
of fermion flavors (N — 2 for polyacetylene). In Eq. (|2.4[) we have used Eq. (|2.2[) to define the GN coupling as 

Agn = NnX TLM = 7 %. • (2.5) 
tori 

III. THE GROSS-NEVEU MODEL IN THE OPTIMIZED PERTURBATION THEORY 

Let us now turn our attention to the implementation of the OPT procedure (20| (for a long, but far from complete 
list of references, please see also Ref. [2l[ and references in there) within the model Lagrangian density given by Eq. 
(|2.4|) . Applying the usual OPT interpolation prescription to the original four- Fermi theory, Eq. (|2.4| . we define the 
interpolated theory 

Cs(ip,ip) = i> (ihlodt - ihvF7id x ) ip — r){l — 8)4>%p + S^^-fivF^ip) 2 , (3.1) 

where rj is an arbitrary mass parameter. It is easy to verify that at S = we have a theory of free fermions, and the 
original theory is recovered for 6=1. Now, by re- introducing the scalar field A, which can be achieved by adding the 
quadratic term (corresponding to a Hubbard-Stratonovich transformation) 

-iS^( A+ ^ hv ^) 2 ' (3 - 2) 

to £s(ip, ip), one obtains the interpolated model corresponding to the original TLM model given by Eq. (|2.3p . 

C s =4> (ih^odt - ihv F jid x ) ip - SAipip - 77(1 - 6)ipip - — ^ N — A 2 . (3.3) 

ZhVFAGN 

Since Eq. (|3 . 3|) is the same model already studied in Ref. [21], so we do not repeat all the details related to the free 
energy density derivation here, where only the main steps and results relevant for our application to the polyacetylene 
will be presented. 

Generally, the OPT method can be implemented as follows. Any physical quantity, is perturbatively computed 
from the interpolated model, up to some finite order-Zc in 6, which is formally used only as a bookkeeping parameter 
and set to the unity at the end of calculation. But in this process any (perturbative) result at order k in the OPT 
remains rj dependent. This arbitrary (a priori) parameter is then fixed by a variational method that then generates 
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nonperturbative results, in the sense that it resums to all orders a certain class of perturbative contributions through 
self-consistent equations. Such optimization method is known as the principle of minimal sensitivity (PMS) and 
amounts to require that be evaluated at the point where it is less sensitive to this parameter. This criterion 
translates into the variational relation 13511 



5=l,f]=rj 







(3.4) 



The optimum value fj that satisfies Eq. (|3.4[) must be a function of the original parameters, including the couplings, 
thus generating "non-perturbative" results. In our case, we are interested in evaluating the optimized free energy at 
finite temperature and density for the scalar field, A, once the fermions have been integrated out. 



A. The Optimized Free Energy Density 

To order-5, Landau's free energy density (or effective potential, in the language of quantum field theories) was 
evaluated in Ref. [21] using functional and diagrammatic techniques. The result is 



■^(Ac, r), T, fx) 



NA 2 



N 





'1 








-(f)] 


2 " 





2(kT) 2 h{r),[i,T) 



TTVpfl 



In 



ln( — ) -7 2 (r/,M,T) 



+ (kT) 2 I 2 3 (r ll ^T) 

where k is the Boltzmann constant and the functions Ii, I 2 and 1% are given respectively by 



(3.5) 



h{r},^,T) = I d.r 



I In 1 + e -V* 2 +W(fcr) 2 -/V(fcT)" + ln \ + e -^Jxi+r?/{kTy+n/{kT) I 



(3.6) 
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h(v, V,T) = / 
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d.r 



^/x 2 + r 1 2 /(kT) 2 



(3.7) 



and 



h (V, M; T) = sinh f A ) / dx , . (3.8) 



\kT 



o cosh (y X 2 +rj 2 /(kT) 2 ^) + cosh [fJ,/(kT)} 



In Eq. (|3 . 5[> . A c is a constant field configuration for the scalar field and M is an arbitrary energy scale introduced 
during the regularization process used to compute the appropriate momentum integrals. In the computation performed 
in Ref. [2l| . the free energy density has been renormalized using the MS scheme for dimensional regularization. We 
also note that Eq. (|3.5p , evaluated at first order in the OPT, already takes into account corrections beyond the 
large- N result. 

By optimizing Eq. (|3.5p through the PMS condition, Eq. (|3.4p . we obtain the optimum value, fj, for the mass 
parameter, which is then re- inserted back in Eq. (|3.5p . allowing us to compute the order parameter A c that minimizes 
the free energy. Using the PMS procedure we then obtain, from Eq. (|3.5p . the general factorized result [2l| 



y(r},iM,T) + rj—y(ri,ii,T) 



(fcT) 2 A GN _ 



0, (3.9) 



where we have defined the function 
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y( v , fi, t) = in (-j - i 2 ( v , t) . (3.10) 

Considering the Xqn/N dependent solution one notices that, when N — > oo in Eq. (|3.9|) , fj — A c and the mean-field 
standard result is exactly reproduced as usual [2l|, [3(|. For finite N, as is our interest here, fj and A c have to be found 
self-consistently by solving the gap equation dT '/dA c — and the PMS equation dT '/drj = (2lj . 



B. The density at finite temperature 



The thermodynamical potential (per volume) is defined as the free energy density at its minimum, f2(T, (i) = 
T(fj, A c , T, fi) and the pressure follows as P(T,fi) — — 0(T, /i). The density is then obtained by the usual relation 
p = dP/dfi. We must also recall that dT jdA c = at A c = A c , due to the gap equation, and that dT '/drj = at 
■q — fj, due to the PMS equation. Then, terms like {dT / dA c )(dA c / dp) and (dT/dr])(dr)/d(J,) do not contribute. One 
then obtains 



p{T,n) = — - 

VFll 



{kTf-I'^^T) +77(77 - A c ) -I^V, l*,T) 



\GN 2 



A 



2tt 2 



V 2 I 2 { V , fx, T)!^, p, T) - ^-{kTfhin, fi, T)I' z {r,, fi, T) 



r;=f;,A c = A c 



(3.11) 



where the primes indicate derivatives with respect to fi. This result will be considered later when we investigate 
thermal effects in y c . 



IV. THE GAP ENERGY AND PARAMETER SET AT T = AND fi = 

In order to perform a numerical analysis we must fix all parameters. This can be done by considering the gap 
energy. In the GN language the order parameter A n is just the TLM gap parameter which, at T = and fi = 0, we 
denote as A . At order-<5 this quantity is given by [2l[ 



where M is an arbitrary (at the moment) renormalization scale to be discussed further below. Eq. (|4.ip explicitly 
includes corrections beyond large-iV, as obtained from our OPT approach. More precisely, taking the mean-field 
approximation, N — * 00 in Eq. (|4.1[) and using the relation A GN = A^ttAtlm, the OPT result exactly recovers the 
mean- field result for TV = 2 [||, 



A MF = M exp[-l/(2A TL M)] , (4.2) 

as one expects [U [3(| • 

Now some remarks concerning the arbitrary energy scale, M, and more generally on the interpretation of Eq. (|4.ip 
in the present polyacetylene context are useful. Usually, in a renormalizable quantum field theory, one can choose 
arbitrary value for M and Aqn will run with the scale appropriately, at a given perturbative order, so that A remains 
scale-invariant as dictated by the renormalization group. For the above gap equation this means that 

-J_ = 1 

Aqn(A^) A G n(Mo) tt V^o 




7 



where Mo is some reference (input) scale 1 . Equations (|4.1[) , or Eq. (|4.3|) , indicates that A(M), or equivalently Ao, 
is the only parameter to be fixed. These equations also show that A(M) — * as In -1 (M/M ) when M — + oo which 
is nothing else than the asymptotic freedom displayed by the GN model. However, in the polymer physics case, the 
interpretation is somewhat different mainly because both Aq and the coupling Aqn are measurable quantities, as we 
shall exploit below. Moreover in contrast to the renormalizable field theory case all quantities here are expected to be 
directly finite, i.e., without need of renormalization due to the explicit high energy cutoff A provided by the 7r-band 
width, i.e., A ~ W. While in our calculation we have used dimensional regularization and renormalization mainly 
for convenience, Eq. (l4.f p should be interpreted as giving the finite, iV-dependent, corrections to the large- N results, 
with the arbitrary scale M (originating from dimensional regularization) to be traded for an explicit cutoff: M = A 
of order A ~ W. Now since M is a parameter from the theory, its precise value is thus a matter of choice to some 
extent, as it does not need to coincide exactly with the experimental parameter W. This implies in particular that 
the scale M can be dealt with in alternative ways as we shall discuss next. 

Consequently, we can consider different possible prescriptions for the basic parameters of the problem, given also 
that some data appear to have non negligible experimental uncertainties. 

• i) First, in the prescription we label (1), Aqn(-W) can be simply set to its phenomenological value given by 



8Na 2 

Agn(m) = W ' (44) 

where we have used the relation 4to = W. As discussed above this implicitly defines a scale M once assuming 
the theoretical prediction of Eq. (|4.f p . 

Regarding the data numerical values, we note that this has been debated for long and different set of values 
appeared in the literature (see e.g. Refs. @,S|). As far as we are aware, it appears however that the present 
widely accepted data values are: K = 21 eV/A 2 , a = 4.1cV/l, 2A = 1.4 - 1.8 eV and W = 4i = 10 eV, 
which are essentially the conventional SSH values [1] except for possible higher values of A [1] , which appears 
as the less accurately determined experimental input. Consequently in our study we shall take this set of input 
but taking the two extreme values of A , that we will call set A and B, respectively for 2A = 1.4 (1.8) eV. 

Comparing thus set A and B appears to us as a very conservative way of taking into account those experimental 
uncertainties, although the higher value of 2Ao ~ 1.8 eV appears to be much more favored in the recent 
literature. Since all relevant physical quantities (such as the critical density) will depend on the cutoff scale, 
M, as already mentioned one possible prescription is to use Eq. (|4.ip to fix the cutoff M value for given 
Xqn(M) = 8Na 2 /(WK) and Ao within accuracy, i.e., for each data set A and B. Eq. (|4.ip shows that due 
to the presence of N dependent constant term 1 — 1/(2 AT), the OPT and the MFA will predict in this way 
different M values even when using the same set of input data parameters. 



• ii) Another possible prescription, that we dub (II), is to set M = W exactly, cutting off the spectrum at an 
energy scale of —W/2 @, [l||. In such case Agn(-^ = W) does not exactly match the experimental value as 
predicted by Eq. (|4.3p . Most previous authors appear to haven chosen this prescription, i.e., changing the 
coupling Aqn value in order that the mean field model best fits the polyacetylene data. We find however equally 
motivated to use the first interpretation since, as already mentioned, the polyacetylene data provides us with a 
rather precise Aqn value, while there is some intrinsic arbitrariness in the precise cutoff scale value (equivalently 
in this case, Ao fixes the energy cutoff scale M within some accuracy) 2 . Actually the two prescription are not 
fundamentally different: in the first one uses the arbitrariness of the cutoff to fit the data Ao and Aqn, while 
in the second one forces the coupling to fit the two scales Ao and W, but this is essentially translating the 
arbitrariness of the scale inside the exponential of Eq. (|4.ip . 



• iii) Finally let us consider yet another possible prescription (or rather interpretation) of the GN 
model/polyacetylene data connexion. It will define our prescription III. Namely, bearing in mind that the 



1 Note that the TV-dependence of the running coupling as dictated by Eq. (I4.3D differs from the standard RG one, which (at one-loop 
order) has a coefficient given by 1 — l/N instead of f — 1/(27V). This difference is a result of OPT which modifies standard one-loop 
order results, in a way expected to improve the latter. 

2 Note also that Eq. (1 4 . . j I) also allows for any other intermediate prescriptions in which M W and \qn(M) ^ 8Na 2 /(W K), but the 
simultaneous equalities for these quantities is excluded because of the 2Aq = f.4 — f .8 eV values. 
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equivalence between the original TLM and continuous GN model was strictly established only at the mean 
field theory level, we may redefine our OPT corrections in the framework of an effective mean field (EMF) GN 
description: more precisely, the OPT-modified gap energy Eq. (|4.1j) can be fitted by the corresponding mean 
field expression Eq. (|4.2[) provided that ones redefines "effective" mean field coupling 3 \* EMF and cutoff W FMF : 



\% MF ee A GN f 1 - 2^) . W emf ^ M { 1 -^) > ( 4 - 5 ) 

together with the identification of these effective parameters to the measured data. 

This freedom of prescriptions, as discriminated above as prescriptions I, 11 and III, actually reflects that neither the 
MFA nor the OPT-improved expression of the gap energy are expected to be exact results. If available, an exact, truly 
non-perturbative calculation of the gap energy would be expected to fit nicely the three independent experimental 
measurements, to (equivalently W), Aq, and Aqn (of course up to limited experimental accuracy). Therefore for 
completeness and comparison purpose, we will consider in the numerical results all these prescriptions together with 
data sets A and B. A summary of the different M and Aqn values for each prescription and data set is given in Table 

in 

Inspection of Table I indeed indicates rather different values of the "bare" coupling Aqn for the three prescriptions, 
which is essentially due to the large uncertainty in Ao between sets A and B. One should not conclude from this that 
our description is lacking prediction. In fact, as we shall see later, the predictions for our main result on the critical 
dopant estimate are not strongly dependent on the coupling values, and will be only slightly different for the three 
cases (provided one uses the same experimental data input). Again, the most important variation will be due to the 
large uncertainty on A . 





LA 


LB 


II.A 


II.B 


III.A 


III.B 


\ OPT 
^GN 


1.28 


1.28 


1.42 


1.55 


1.89 


2.07 


\MFA 
^GN 


1.28 


1.28 


1.18 


1.30 


1.18 


1.30 


M OPT (eV) 


13.85 


17.80 


10 


10 


7.5 


7.5 


M MFA (eV) 


8.15 


10.47 


10 


10 


10 


10 



TABLE I: The OPT and MFA values for Agn and M obtained, from the gap equation, for parameter sets A (2Ao = 1.4 eV) 
and B (2Ao = 1.8 eV), using prescriptions I— III for fixing the relevant parameters of the model. The common values for both 
cases are: K = 21 eV/A 2 , a = 4.1 eV/A and W = 4t = 10 eV. 



V. PHASE DIAGRAMS 



Having set the parameters for different prescriptions we can investigate the phase diagrams for the theory. Let us 
start by locating the second order and first order transition lines in the T — fj, plane. This is shown in Fig. [T] for 
the choice of prescription IB. It shows the appearance of a tricritical point around kT/M ~ 0.012 and /i/M ~ 0.025. 
Those numbers would slightly change for the other prescriptions, with the overall behavior qualitatively very similar. 
Although the appearance of a tricritical point is an interesting issue when considering the GN as a toy model for QCD, 
it has no practical implications for the polyacetylene, in which case one is concerned with temperatures lower than 
about Td ~ 400 K , above which the polyacetylene is unstable and decomposes when heated (instead of melting) [33] , 
that is kTjIM w 0.0020 for M « 17.80 eV (prescription IB). This type of phase diagram has been extensively studied 
in Ref. [211 ], where analytical expressions for T c , p c , as well as useful relations among tricritical points relations can 
be found. 

As emphasized in the introduction, one of our goals here regards the evaluation of the critical concentration, y c , 
as a function of the temperature. With that aim one benefits from analyzing the phase diagram in the T — p plane 
since y c (T) is directly proportional to p c {T). This is shown in Fig. [5] The dot in Fig. [5] indicates the tricritical 
point above which the transition is of the second kind. The mixed (semi conductor) region is associated to the 



3 Note that the meaning of effective coupling here is purely phenomenological, as obtained from a fit, and thus unrelated to the usual 
effective coupling Aqn(M) as above discussed having the RG behavior. 
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kT/M om 



0.015 



0.02 




(CSB) ^ 

\ 

\ 



0.005 



\ 











0.005 0.01 0.015 



0.02 0.025 



FIG. 1: The OPT phase diagram, in the kT/M-p/M plane, for TV = 2 and A = 1.28 (prescriptions IA and IB in Table [Q. The 
continuous line represents the second order transition whereas the dashed line represents the first order transition and the dot 
indicates the tricritical point, which occurs ate kT/M ~ 0.012 and p/M ~ 0.025. 

first order phase transition. At T — the critical dopant density is approximately p c (0) — 0.016 M/(vpK). Figure 
[5] shows the situation where the first order transition line, which appears in the T — p plane, splits into two lines 
limiting a coexistence, mixed (semiconducting) region. We note that there are indeed experimental indications of 
a mixed phase for polyacetylene for concentrations below the critical one [38|. It is also interesting to note, from 
the same figure, that when one evaluates y c at T = using the GN model [17l. HH the only observed transition is 
from the semiconducting phase to the metallic one. However, even at room temperature (roughly kT/M w 0.0015) 
our figure displays another transition from the (unsymmetric) insulator phase to the (mixed) semiconductor phase 
which happens at a rather very small density, of the order vf^p/M ~ 10 -5 . For this transition the critical density 
increases with the temperature. On the other hand, the critical density when going from the (mixed) semiconducting 
phase to the (symmetric) metallic one seems to slightly decrease for low values of T. Figure also shows that above 
p c (0), computed in the next section, the material is a conductor at any temperature, provided that this temperature 
is smaller than the degradation temperature, Td- In the next section we shall devote especial attention to this issue, 
since the literature does not seem to indicate any previous studies of the influence of thermal effects in y c within the 
models considered here. 












0.005 



5 0.01 

v F h p/M 



0.015 



FIG. 2: The OPT phase diagram in the kT/M - v F hp/M plane for TV = 2 and A = 1.28 (prescriptions IA and IB in Table 
The insulator region is associated with the unsymmetric (dimerized) phase, while the metallic region is associated with the 
symmetric (undimerized) phase. 
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VI. CRITICAL DOPANT DENSITY 

In this section we evaluate the critical dopant concentration, y c , which, with periodic boundary conditions in the 
polyacetylene chain, is given simply as y c = ap c . In the next subsection we consider the case T = performing a 
numerical comparison between the MFA and the OPT using the different sets of parameters presented earlier. Next, 
we will consider how the temperature affects y c . 

A. The zero temperature case 

In this section we will evaluate the dopant critical density in polyacetylene within the OPT approach, neglecting 
eventual temperature effects which will be considered in the next subsection. Let us start by taking the limit T —> 
in the free energy density, Eq. (|3.5[) . The various functions defined by Eqs. (|3.6p , (|3.7[) and (|3.8[) . in the T — » limit 
become 



Urn (kT) 2 h( V ,p,T) = -±6(p-r,) 



2, f M + vV - V 2 I — 

r\ In - py p z - V 



(6.1) 



lim h (v, Mi T) = 9(p - r,) In ( M + V ' ^ ^ ) , (6.2) 



r— o 



»7 



UmfcT7 3 (»7,/i,T) = sgn^)^ - 77) vV - ?7 2 • (6.3) 

Using Eqs. (|6.1[) . (|6.2p and (|6.3p in Landau's free energy density, Eq. (|3.5[) . we notice that it can be divided into two 
cases: i) p < r\ and ii) p > 77. At zero temperature the critical chemical potential, ^ c (0), is defined as the one which 
produces the same pressure for both, A c = Ao 7^ and A c = 0. This quantity, which has been evaluated in Ref. [2l| . 
is given by 





»H-S^rjMi-^) . (6 .4) 

As expected, there are two values of p corresponding to /z c (0). The first is simply p — 0, corresponding to the minimum 
of the free energy density that occurs at A c = Ao 7^ for the case p < rj (corresponding to the T = p, = Q situation). 
In the second case (p > 77), the minimum of the free energy density occurs at the origin, A c = 0. In this case the 
PMS relation, Eq. (|3.9p . implies that fj = and one gets, from Eq. (|3.11|) and after simple algebra, the result 

p c (0)N f A GN \ 
PM = —r— 1 - 7T-T7 , (6-5) 



itTivf \ 2ttN _ 

where the multiplicative factor in the RHS of the above equation, p c (0)N/ (ttTivf), is just the MFA result (with 
N = 2). In Eq. (|6.5[) the term Agn/(27t./V) gives the first order OPT finite N corrections to /O c (0). One can now 
insert Eq. (|6 .4[) into Eq. (|6.5[) . Using also the expression for Ao, Eq. (|4.ip . we obtain an analytical expression for the 
critical density that includes finite N corrections, 

„„ ( ) = aM 0) = ^( 1 -^LV / 7l-^, (6-6) 



nW V 2ttA/ V 2N 



where we have used the defining relation cl/(Tlvf) = 2/W. Note that the explicit (overall) scale M dependence 
disappears from Eq. (|6.6p . where only the physical parameters W/Aq set the overall scale. However, there is an 
implicit scale dependence in the OPT case via the dependence on Aqn(M) within the factor [1 — Agn/(27tA)] 1 ^ 2 4 . 



4 Note, however, that since [1 - A G N/(27r-/V)] 1/2 [1 - 1/(2JV)] — > 1, as N — > 00, the MFA results only depend on the ratio 2A /W 
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LA 


LB 


II.A 


II.B 


III.A 


III.B 


y? PT 


0.0448 


0.0576 


0.0450 


0.0569 


0.058 


0.074 


yf FA 


0.0630 


0.081 


0.0630 


0.081 


0.0630 


0.081 


t/a 


7.143 


5.555 


7.143 


5.555 


7.143 


5.555 



TABLE II: The critical dopant concentration y c obtained with the OPT and the MFA obtained from each of the three different 
prescriptions and two data sets at T = 0. For reference we also show the relevant values of £/a in each case. 

Concerning our third prescription, the "effective mean field" (EMF) prescription III defined in Sec. IV, where the 
OPT corrections are reinterpreted differently as redefining effective coupling and scale, the corresponding expression 
of the critical dopant concentration is straightforward to derive using definitions Eq. (|4.5p and reads: 



,EMF 



(0) 



V2NAo 



nW, 



EMF 



^*EMF 



1/2 



(6.7) 



where it is again understood that corresponding set A or B data values should be used now for X* EMF and W FMF . 

We are now in position to make predictions concerning the observable y c (0) using, for completeness, the two different 
sets of data parameter and our three different theoretical prescriptions. The comparison between the OPT and MFA 
results is shown in Table II. 

As one can see from Table II, the results depend rather substantially on the experimental data set choice, due 
essentially to the (linear) dependence on Ao. The most quoted experimental value of y c is y c ~ 0.06, although its 
precise value is not very accurately determined. Typically, looking e.g. at the data from Refs. (?l [3^|. one may infer 
non negligible uncertainties on the exact transition value, of about Ay c ~ 0.01. Also, slightly higher values of y c have 
been reported in other studies [40l |. In contrast, the two different prescriptions I and II regarding the scale dependence 
only affect mildly the y c values in the OPT case, i.e., the arbitrary scale dependence, which appears only indirectly 
within the factor [1 — Xq-^/(2t:N)], which appears in our previous expressions like Eq. (|6.6|) . remains moderate. 

Inspection of Table II indicates that the OPT performs better for smaller values of £/a. We have carried out 
numerical simulations that indicate that, in fact, the OPT and MFA predict similar deviations from the experimental 
value, y c — 0.06, for £/a ~ 6.4. For £/a < 6.4, however, the OPT predictions are better than the MFA ones and the 
situations gets reversed for £/a > 6.4 as Table II shows. To illustrate this point we present Fig. [SJ where y c (0) is 
plo tted as a function of £/<x, using prescription IB (Tabic IJ). This figure shows that for the value £/a = 6, which Ref. 
[19| refers to, the OPT (y c = 0.0535) performs better than the MFA (y c = 0.0750). We also remark that the same 
pattern is obtained when one considers others sets of values and prescriptions. 



y c (0) 



0.08 
0.07 
0.06 
0.05 



5.5 



6.5 



7.5 



FIG. 3: The OPT (continuous line) and MFA (dashed line) predictions for j/ c (0) as a function of £/a, using prescription IB 
(Table |IJ, in the relevant range 5.5 < £/a < 7.5. The experimental value is y c (0) — 0.06. 



Finally, as far the different prescriptions are concerned, we note that the ones that make use of the band gap energy 
as 1.8 eV (used in our data sets B in Table I), produces results with much better agreement with the experimental 
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data for the OPT than the MFA, as one can check from the results in Table II. In fact, this value for the band gap 
energy appears in the literature as a more satisfactory value for polyacetylene Q . 

B. The Finite Temperature Case 

To obtain the temperature dependence of the critical dopant concentration one can proceed numerically considering 
y c (T) = ap c (T). First, as already emphasized, one should note that in practice there is an upper temperature of about 
Td ~ AQQK above which our simple models break down since the polymer undergoes other phase transitions. Figured] 
shows that from the absolute zero temperature to the upper temperature our prediction to the decrease in the critical 
dopant concentration is only about 1%, while it is less than 0.5% from room temperature (about 300-ftT), where most 
of the experiments are done, to the upper temperature. This shows that in practice, at least with the type of models 
considered here, one may safely evaluate y c at T = as it has been done in Refs. [UGH- 



0.0576 
0.0575 
0.0574 

y c 

0.0573 
0.0572 



100 200 300 400 500 

T(K) 

FIG. 4: The critical dopant concentration, y c (T), temperature dependence for parameters set B and procedure I. The degra- 
dation temperature is about Td ~ 400K. 



VII. CONCLUSIONS 

Using the OPT we have reviewed the evaluation of Landau's free energy for the 1+1 dimensional massless Gross- 
Neveu model at finite temperature and chemical potential as performed in Ref. [2l[ . Then, relating this model to the 
TLM continuous model for polyacetylene, we have computed the critical dopant concentration, j/ c , for the transition 
to the metallic phase. The (divergent) free energy density has been rendered finite by using a renormalization 
procedure which is standard in quantum field theories (MS scheme with dimensional regularization). An arbitrary 
energy scale, M, introduced during the formal regularization process was fixed by using polyacetylene experimental 
inputs. Regarding the matching of the theory parameters to the experimental data values, we have provided a critical 
discussion of the possible different prescriptions to relate them giving, at the same time, the expected results for each 
of the prescriptions used. 

The OPT formalism allows for the inclusion of finite N corrections already at the first non trivial order which should 
improve the usual MFA results since, for this particular case, N = 2. To illustrate the possible phase transitions 
allowed by the GN model we have obtained phase diagrams in the T — p as well as in the T — p planes. Then, we have 
obtained a neat analytical expression for y c = ap c at T = which contains explicit 1/N corrections. Our results show 
that when one uses up-to-date parameter values, the OPT results improve over the MFA as expected and are in good 
agreement with the experimental result. Another result of the present work regards the study of possible thermal 
effects in y c . Our analysis has been performed in a numerical fashion, showing that, for a realistic temperature range, 
< T < 400 K, thermal effects induce a negligible decrease of y c when going from the (mixed) semiconducting phase 
to the (symmetric) metallic phase. 

Regarding the metal-insulator transition, one must recall the importance of the transport property and the local- 
ization problem [4ll. |42|. As far as polymers are concerned, it is known that their transport properties result from 
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the mechanism of hopping [421] . which leads the conductivity to increase with the temperature till some maximum 
value. At the same time, disorder in the system can result in the localization of states and, if it is too strong, it can 
lead to an insulator behavior. The metal-insulator transition in real systems is then a consequence of the interplay of 
the amount of disorder, doping and the thermal activation process. On the theoretical side, an interesting possibility 
of extending the field theory application method we used to study the metal-insulator transition in polyacetylene, 
would be the calculation of the conductivity a and the determination of its T, fi dependence. This would then allow a 
closer comparison with the experimental measures on this quantity. Basically, from response theory, the conductivity 
can be computed from a Green-Kubo formula, which entails a calculation of specific correlation functions and higher 
loop Feynman diagram contributions in the GN model we used. Although the calculation of a in this context has 
been already considered within other approaches and approximations in the literature, we hope to pursue a detailed 
calculation of <r(T, /z, N) including finite N corrections, which we believe has not been done up to now. However, this 
is a non trivial calculation which is well beyond the scope of the present paper. We intend to address this issue in 
the future. Another interesting possible extension of the present work would be to use the OPT and the methods 
developed in Ref. [34| to consider the massive GN model, which can be related to the czs-polyacetylene. 
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